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ABSTRACT 

It has been pointed out that the Galactic satellites all have a common mass around 10 7 M Q within 300 pc (M0.3), 
while they span almost four order of magnitudes in luminosity (Mateo et al. 1993, Strigari et al. 2008). It is 
argued that this may reflect a specific scale for galaxy formation or a scale for dark matter clustering. Here we 
use numerical simulations coupled with a semi-analytic model for galaxy formation, to predict the central mass 
and luminosity of galactic satellites in the ACDM model. We show that this common mass scale can be explained 
within the Cold Dark Matter scenario when the physics of galaxy formation is taken into account. The narrow 
range of Mo .3 comes from the narrow distribution of circular velocities at time of accretion (peaking around 20 
km/s) for satellites able to form stars and the not tight correlation between halo concentration and circular velocity. 
The wide range of satellite luminosities is due to a combination of the mass at time of accretion and the broad 
distribution of accretion redshifts for a given mass. This causes the satellites baryonic content to be suppressed 
by photo-ionization to very different extents. Our results favor the argument that the common mass Mo .3 reflects 
a specific scale (circular velocity ~ 20 km/s ) for star formation. 

Subject headings: cosmology: theory — galaxies: dwarf — hydrodynamics — methods: numerical 



1. INTRODUCTION 

Galaxies are thought to form out of gas which cools and col- 
lapses to the center of dark matter haloes (White & Rees 1978). 
A correlation between stellar mass and host halo mass is thus 
expected. Observational data, such as from the SDSS, have 
shown a good correlation between the stellar mass of central 
galaxies and the halo mass of galaxy groups (e.g. Yang et al. 
2007), but such a correlation is not mantained in satellite galax- 
ies (e.g., Gao et al. 2004, Conroy et al. 2007). It is also not 
clear if such a relation persists to all mass scales. For exam- 
ple, in galaxy clusters, the stellar mass of the central galaxy is 
not a strong function of the halo mass, because gas cooling and 
star formation are regulated by the physical processes operating 
in clusters, such as AGN feedback, which gives rise to a steep 
luminosity function at the bright end (e.g., Kang et al. 2006). 

Mateo et al. (1993) and more recently, Strigari et al. (2008, 
hereafter S08), claimed that all the satellites with luminosity 
between 10 3 and 10 7 Lq to have a common mass of ~ 10 7 M Q 
within a radius of 300pc (Mo 3), and to be dark matter domi- 
nated within this region. This narrow range for central masses 
compared to the wide range of luminosities has raised the is- 
sue if this common mass may reflect a specific scale for galaxy 
formation or a scale for dark matter clustering. 

In this letter, we use a series of high-resolution simulations 
of MW type haloes combined with a semi-analytical model of 
galaxy formation to show how this flat relation between M0.3 
and luminosity arises naturally within our standard cosmolog- 
ical framework combined with simple modelling of the galaxy 
formation process. 



2. SIMULATIONS 

Nbody simulations have been carried out using PKDGRAV, a 
treecode written by Joachim Stadel and Thomas Quinn (Stadel 
2001). The cosmological parameters are set in agreement with 
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the WMAP mission first year results (WMAP1 : Spergel et al. 2003) 
as follows: O\=0.732, ft m =0.268, 14=0.044, h = 0.71, n = 1.0 
and a% = 0.9. We have selected three candidate haloes with a 
mass similar to the mass of our Galaxy (M - 10 I2 M Q ) from 
an existing low resolution dark matter simulation (300 3 parti- 
cles within 90 Mpc) and re-simulated them at higher resolu- 
tion using the volume renormalization technique (Katz & White 
1993). Our high resolution haloes all have a quiet merging his- 
tory with no major merger after z = 2, thus are likely to host 
a disk galaxy at the present time (as confirmed by our semi- 
analytic models, see section The high resolution run is 12 3 
times better resolved than the low resolution one: the dark mat- 
ter particle mass is =4.16 x 10 5 h~ l M Q , where each dark 
matter particle has a spline gravitation (comoving) softening of 
355 h~ l pc. Properties of single halos are listed in TableQ] 

For the purpose of constructing accurate merger trees for 
each simulated galaxy we analyzed 53 outputs between z = 20 
and z = 0. For each snapshot we looked for all virialized haloes 
within the high resolution region using a spherical overden- 
sity algorithm (see Maccio et al. 2007 for more details on our 
halo finding procedure). We included in the halo catalogue all 
haloes with more than 100 particles (> 4 x 10 7 /r'M Q ). For 
the merger tree construction we started marking all the parti- 
cles within the 1 .5 times the virial radius of a given galaxy at 
z = and we tracked them back to the previous output time. 
We then make a list of all haloes that at earlier output time 
containing marked particles, recording the number of marked 
particles. In addition we record the number of particles that 
are not in any halo in the previous output time and we consider 
them as smoothly accreted. We used the two criteria suggested 
in Wechsler et al. (2002) for halo 1 at one output time to be 
labeled a "progenitor" of halo 2 at the subsequent output time. 
In our language halo 2 will then be labeled as a "descendant" 
of halo 1 if i) more than 50% of the particles in halo 1 end up 
in halo 2 or if ii) more than 75% of halo 1 particles that end 
up in any halo at time step 2 do end up in halo 2 (this second 
criterion is mainly relevant during major mergers). A halo can 
have only one descendant but there is no limit to the number of 
progenitors. On average there are 20,000 progenitors for each 
high resolution DM halo. 
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3. SEMI ANALYTICAL MODELS 

We calculate the luminosity of the satellites using the semi- 
analytical model for galaxy formation of Kang et al. (2005), 
and we refer the reader to this paper for more details. Basi- 
cally we follow the formation history of dark matter haloes and 
graft the physical models for galaxy formation onto the merger 
trees. Physical processes governing galaxy formation include: 
gas cooling from radiative hot gas, star formation and super- 
nova feedback. The original model has been recently updated 
(Kang 2008) to include photo-ionization to suppress baryon ac- 
cretion in low mass haloes (Kravtsov et al. 2004), along with 
a new fitting formula to describe dynamical friction time-scales 
for satellite galaxies (Jiang et al. 2008). Kang (2008) has 
shown that this model is successful in reproducing the Milky 
Way satellite luminosity function and other properties (Maccio 
et al. in preparation). 

There are two important factors which regulate star forma- 
tion efficiency in low mass dark matter haloes. Firstly, the cool- 
ing rate for hot gas in haloes with a virial temperature T v / r < 
10 4 K is quite low due to the inefficiency of H2 cooling, thus we 
shut off gas cooling in such haloes. This effect is responsible for 
the absence of visible satellites with y c i lc < 16.6 km/s as shown 
in Figure [3] Secondly, hot gas accretion is suppressed by the 
cosmic photo-ionization background in low mass haloes; the 
typical filtering-mass, at which haloes lose half of their bary- 
onic content, has a strong redshift dependence and it increases 
from 1O 7 M at z ~ 10 to 3 x 1O 1O M at z = (Kravtsov et al. 
2004). This implies that a satellite galaxy with a halo mass 
around 10 9 M Q at infall, has its baryon content suppressed to 
different extents depending on its accretion redshift. We will 
see in Section|4]that this gives rise to the wide luminosity range 
for satellites, despite the narrow range spanned by circular ve- 
locities at the accretion time. We set the reionization redshift to 
be z r = 7 but as shown by Kravtsov et al. (2004) the results of 
photo-ionization are almost insensitive to the actual value of z r - 

4. RESULTS 

The relation between the central mass of a satellite and its 
circular velocity can be predicted on a theoretical basis, using 
the well studied relation between mass and concentration. N- 
body simulations have shown that the spherically averaged den- 
sity profiles of dark matter haloes can be well described by a 
two parameter analytic profile (Navarro, Frenk & White 1997, 
NFW hereafter): 

^2 = *1 (1) 

Pcrit (r/r s )(l + r/r s ) 2 ' 

where p CI j t is the critical density of the universe, 8 C is the char- 
acteristic overdensity of the halo, and r s is the radius where the 
logarithmic slope of the halo density profile dlnp/dlnr = -2. 
A more useful parametrization is in terms of the virial mass, 



Table 1 
Galaxies Parameters 



Halo 


Mass 


Npart 


Rvir 


MaxVcirc 




(IO^Mq) 




(kpc/h) 


(km/s) 


GO 


0.88 


2115385 


197 


179.6 


Gl 


1.22 


2931295 


219 


187.6 


G2 


1.30 


3123511 


250 


203.1 



M v ir, and concentration parameter, c = R v - u /r s . The virial mass 
and radius are related by M v ; r = A v i r /5 cl -it(47r/3)7?^ ir , where A v jr 
is the density contrast of the halo. Then the mass within 300 pc 
is simply defined by: 

M 03 = / 2T7t 2 p(t)dt = 4-Trp crit S c rl(log(l+x)-- ) (2) 

Jo ' l+x 

where x = ro .3 / r s and the value of M0.3 only depends on the two 
parameters defining the density profile: r s (or the concentration 
c) and S c . Several models have been proposed to link mass and 
concentration of dark matter haloes (Bullock et al. 2001, Eke 
et al. 2001, Maccio et al. 2007), recently Maccio et al. 2008 
have proposed a new toy model (based on a modification of 
original model of Bullock et al. 2001) to predict the concentra- 
tion mass relation in dark matter haloes. Using this model it is 
then possible to predict the value of Mo. 3 as a function of the 
circular velocity of the halo. 

Figure Q] shows the values for the mass within 300 pc in the 
WMAP5 cosmology (Komatsu et al. 2008) as a function of the 
circular velocity of the dark matter halo. Dashed lines represent 
the one a scatter as expected from Nbody simulations (here we 
use a constant value of 0.1 dex Maccio et al. 2008). The dotted 
horizontal lines show the mass range observed for Milky Way 
satellites (S08). The range of circular velocities that is in agree- 
ment with observational data (shaded gray region in the figure) 
is fairly large and covers almost the entire range of plausible 
values for V c i lc . 

There are two possible issues with this kind of approach: 
firstly, we are extrapolating to lower masses a relation (c/M) 
that has been tested on N-body simulation only down to a mass 
of ~ 10 10 Mq, well above the expected masses for Milky Way 




V circ (z = 0) [km/s] 

FIG. 1. — Theoretical prediction of the mass within 300 pc versus circular 
velocity using the toy model of Maccio et al. 2008 for the mass concentration 
relation: solid line and dashed lines represent the median and two <r scatter 
respectively for the WMAP5 model (the dotted line shows the median for the 
WMAP1 model). The dotted (red) horizontal lines show the mass range ob- 
served Milky Way satellites (S08). The shaded region gives a visual impression 
of circular velocities for which the values of Mq .3 consistent with observations 
are expected. 
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satellites. Secondly, satellites we see today can be a biased sam- 
ple of the overall dark matter halo population since they are the 
surviving population which possibly form prior to reionisation 
(Moore et al. 2006). For these reasons a full numerical inspec- 
tion of origin of the narrow range for M0.3 is needed. 

In order to compare numerical results to observations, for 
each simulated satellite we need both its luminosity L and its 
inner mass M0.3. The first quantity is a direct outcome of our 
semi-analytical model, whilst to compute the second quatity we 
proceeded in the following way: at time of accretion of each 
satellite we compute the density profile directly from its particle 
distribution in the N-body simulation. The resulting numerical 
density profile is then fitted with an NFW profile (Eq.[T|i; dur- 
ing the fitting procedure we treat both r s and S c as free param- 
eters. Their values, and associated uncertainties, are obtained 
via a \ 2 minimization procedure (see Maccio et al. 2008 for 
more details). We are only interested in dark matter haloes that 
host a galaxy according to our semi-analytic model; given that 
gas cooling is allowed only in haloes with M > 10 8 /i _l M Q (i.e. 
T V ir > 10 4 K) this implies that, on average, we have more than 
1,000 particles per halo at time of accretion, which is sufficient 
to obtain a robust estimation of the density profile parameters 
(Maccio et al. 2007). Under the assumption that the density pro- 
file within 300 pc does not evolve from the time of accretion to 
Z = we can compute, for each satellite, the present value of 
M0.3 using equation|2] 

The upper panel in Figure [2] shows the results for the rela- 
tion between the mass within 300 pc and luminosity as obtained 
in our numerical model (red dots) versus the observational re- 
sults (black dots with error bars). Here we plot results only for 
simulated satellites that satisfy the detection threshold of the 
SDSS as determined by Koposov et al. 2007. This means that 
the satellite luminosity and distance have to satisfy the follow- 
ing relation log(/?/kpc) < 1.04-0.228 M v . The mean and the 
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FIG. 2. — Mass within 300 pc versus luminosity. Red dots show results from 
our numerical model, black points with error bars are the observational results 
from S08. Upper panel: no correction for the concentration related density 
evolution. Lower panel: correction included (see text). 



scatter of observational data are both well reproduced by our 
numerical results up to a luminosity of L = 2 x 10 6 Lo, after this 
point simulations seem to suggest an increase with luminos- 
ity of M0.3, which is not present in the data (even if only three 
satellite galaxies have a luminosity greater than 10 6 L Q ). 

These results are obtained under the assumption of no evo- 
lution for the parameters defining the density profile (r s and<5 c ). 
This assumption is motivated by the detailed numerical study 
carried out by Kazantzidis et al. 2004 (K04, hereafter, see also 
the recent results by and Penarrubia et al. 2008). K04 have 
shown that the inner density profile is extremely robust and 
that it is unmodified by tidal forces even after tidal stripping re- 
moves a large fraction of the initial mass. They have also shown 
that the degree of modification (if any) of the density profile de- 
pends on the initial (i.e. before infall) concentration of the satel- 
lite dark matter halo. While highly concentrated haloes (with 
c > 15) are able to keep the profile unchanged even after several 
orbital periods, less concentrated haloes (c « 9) slightly modify 
their profile mainly by reducing the overall normalization (S c ) 
by approximately a factor 2 (see also Mayer et al. 2006). 

In order to take into account this expected modification of the 
density profile in low concentration haloes, we manually reduce 
the (5 e parameter by a factor 2 in all haloes with concentration 
less than 10 at the moment of infalling, In doing this correc- 
tion we used the value of the concetration extrapolated to z=0, 
in other words we multiplied the value of the concetration at 
Zacc by (1 +Zacc) in order to take into account the redshift evolu- 
tion of R wu . Results are shown in lower panel of Figure [2] As 
expected this modification mainly applies to high luminosity 
haloes, since they were the most massive ones at time of infall, 
and thus likely to be less concentrated (the average concentra- 
tion of our haloes is 16.3 ±7.1, and 74% have c > 10). When 
a possible modification of the density profile for low concen- 
tration haloes is taken into account the up turn in the numerical 
M0.3/L relation at high luminosities almost vanishes and numer- 
ical results are now in better agreement with the observational 
data. There are still some haloes with Mo. 3 > 3 x 10 7 around 
L = 10 6 , these haloes formed at high redshift, and thus happen, 
by chance, to have a high c (i.e. they are not affected by our 
correction) and a large L; besides of that they do not present 
any other peculiar behavior. 

The presence of a baryonic component inside the dark matter 
halo can by itself modify the density profile of the halo, due to 
the adiabatic compression process (e.g. Blumenthal et al. 1986). 
In our case we expect this effect to be negligible given that the 
baryon fraction of our haloes is much lower then the universal 
one and satellites have been observed to be dark matter domi- 
nated even in their central regions (S08). 

Now we turn to understand the origin of the relation between 
M0.3 and L found by S08. In upper panel of Figure [3] we show 
the distribution of V C!rc at the time of accretion for visible satel- 
lites (same sample used in Figure 0. The distribution peaks 
around Vd lc = 20 km/s and then declines sharply towards higher 
values of the circular velocity. As discussed in Section [3] the 
sharp cutoff below V^c ~ 20 km/s comes from the shut off of 
cooling in haloes with virial temperature below 10 4 K. Combin- 
ing this narrow distribution of circular velocity (between 20-40 
km/s) with the theoretical expectations shown in Figure Q] and 
assuming that Mo .3 does not evolve after accretion, it is not sur- 
prising that all Milky Way satellites (observed and simulated) 
have a inner mass within 300 pc always around 10 7 M Q . It is 
then interesting to ask why these satellites span a wide range 
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of luminosity. In the lower panel of Figure [3] we show the 
dark matter mass of satellites at their time of accretion, and the 
dashed line shows the evolution of the filter mass, defined as the 
mass of a halo in which half of its baryons have been lost due 
to photo-ionization. We can see that most satellites have mass 
lower than the filter mass at accretion, and at a given accretion 
mass, there is wide range of accretion redshifts, which give rise 
to different amounts of baryon suppression. The spread in Lu- 
minosity orginates then from the range of halo masses at time 
of accretion combined with the large scatter (at a given mass) 
in the accretion time. 

Recently Okamoto, Gao & Theuns (2008) have shown that 
the actual values of the filter mass might be smaller than what 
suggested by Kravtsov et al. (2004). The effect of a lower filter- 
ing mass will be to increase the satellite luminosity for a given 
y c i rc ; this will push all points in Figure [2] towards the right; but 
will not alter the flat relation between 1O 4 L0 and 1O 7 L0. 

Whilst this paper was ready for submission a similar study 
was posted on the arXiv (Li et al. 2008). These authors used a 
similar approach (an Nbody simulation combined with a semi- 
analytic model for galaxy formation) to study the relation be- 
tween central mass and luminosity in (satellite) dwarf galax- 
ies. The main difference with respect to this work is that Li 
et al. (2008) computed the inner mass of satellites directly from 
the Nbody simulation. Given the spatial resolution of their sim- 
ulation (0.18 kpc) they presented results for Mo.6 (mass within 
600 kpc) which turned out to be in excellent agreement with 
observational data of Strigari et al. (2007). This study is com- 
plementary to ours. 

5. SUMMARY AND CONCLUSIONS 

The observational evidence that all Milky Way satellites have 
a common mass of about 10 7 M Q within their central 300 par- 
sees has raised issues about the possible existence of a new 
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FIG. 3. — Upper panel: distribution of the circular velocities at the time of 
accretion for visible satellites. Lower panel: accretion redshift vs satellite halo 
mass at accretion time. The (red) dashed line shows the redshift evolution of 
the filtering mass for UV photoionization. 



scale in galaxy formation or a characteristic scale for the clus- 
tering of dark matter. In this Letter, by using high resolution 
numerical simulations combined with a semi analytic model for 
galaxy formation, we show that this common mass scale can be 
easily explained within the current (A)CDM model for structure 
formation. 

The observational data on the L/M0.3 relation can be success- 
fully reproduced in numerical simulations, up to a luminosity 
of 10 6 Lq under the assumption that the parameters describing 
the density profile of satellite galaxies do not evolve signifi- 
cantly after these galaxies have been accreted into the main 
halo. When a plausible (small) modification of such param- 
eters, for low concentration haloes, is taken into account, the 
agreement with observational data extends to the whole lumi- 
nosity range. 

According to our numerical modeling this common mass scale, 
M0.3 ~ 10 7 M Q , originates from the narrow range of circular ve- 
locities (20-40 km/s) spanned by visible satellites at the time of 
accretion (with a lower limit of 17 km/s set by the inefficiency 
of H2 cooling in halo with virial temperature below 10 4 K). On 
the other hand the wide range of luminosities comes from the 
range of halo masses at time of accretion combined with the 
large scatter in the accretion time for an halo with a given mass, 
which has a strong impact on the ability of photo reionization of 
reducing the baryon content of satellites, and thus determining 
their luminosity. 

Our results show that the observed flat relation between satel- 
lite inner mass and luminosity can be easily explained within 
the CDM framework, without modifying the nature of dark 
matter particles. 

The authors thank L. Strigari for making his results publicly 
available. Frank van den Bosch and Sergey Koposov are grate- 
fully acknowledged for interesting discussions and useful com- 
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Planck-Institut fur Astronomie at the Rechenzentrum in Garch- 
ing. 
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